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ABSTRACT 



psj . We present ultraviolet, optical, and infrared photometry of the afterglow of the X- 



ray flash XRF 050416A taken between approximately 100 seconds and 36 days after 
the burst. We find an intrinsic spectral slope between 1930 A and 22 200 A of /? = 
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— 1.14 lb 0.20 and a decay rate of a = —0.86 it 0.15. There is no evidence for a change 
in the decay rate between approximately 0.7 and 4.7 days after the burst. Our data 
imphes that there is no spectral break between the optical and X-ray bands between 0.7 
and 4.7 days after the burst, and is consistent with the cooling break being redward of 
the Kg band (22 200 A) at 0.7 days. The combined ultraviolet /optical/infrared spectral 
energy distribution shows no evidence for a significant amount of extinction in the host 
galaxy along the line of sight to XRF 050416A. Our data suggest that the extragalactic 
extinction along the line of sight to the burst is only approximately Ay = 0.2 mag, which 
is significantly less than the extinction expected from the hydrogen column density 
inferred from X-ray observations of XRF 050416A assuming a dust-to-gas ratio similar 
to what is found for the Milky Way. The observed extinction, however, is consistent with 
the dust-to-gas ratio seen in the Small Magellanic Cloud. We suggest that XRF 050416A 
may have a two-component jet similar to what has been proposed for GRB 030329. If 
this is the case the lack of an observed jet break between 0.7 and 42 days is an illusion 
due to emission from the wide jet dominating the afterglow after approximately 1.5 
days. 

Subject headings: gamma rays: bursts 



1. Introduction 

The Swift satellite (Gehrels et al. 2004) is a multi-instrument observatory designed to detect 
and rapidly localize gamma-ray bursts (GRBs). The observatory contains three instruments. The 
Burst Alert Telescope (BAT; Barthelmy et al. 2005) is used to identify GRBs and localize them 
to < 3' in the energy range 15-150 keV. Once a burst has been detected Swift slews to point 
its two narrow-field instruments, the X-Ray Telescope (XRT; Burrows et al. 2005), and the Ul- 
traViolet/Optical Telescope (UVOT; Roming et al. 2005) at the burst. The XRT obtains X-ray 
localizations to < 5" in the energy range 0.2-10 keV. The UVOT simultaneously obtains localiza- 
tions to ~ 0'.'5 then cycles through a set of optical and ultraviolet filters. 

The X-ray flash (XRF) XRF 050416A was detected in the constellation Coma Berenices by the 
BAT at 11:04:44.5 UTC (Sakamoto et al. 2005a). The gamma-ray light curve showed a slowly rising 
peak followed by several smaller peaks. The burst had a Tqq duration of only 2.4ib0.2 s in the 15-150 
keV band. Sakamoto et al. (2006) find a peak energy of Ep = 15.6+1^ keV, making XRF 050416A 
a soft burst, and they classify XRF 050416A as an XRF. The fluence was (3.5 ± 0.3) x lO"'^ erg 
cm~^ in the 15-150 keV band (Sakamoto et al. 2006). The fluence in the soft 15-25 keV band 
((1.7 ± 0.2) X lO"'^ erg cm'^) was larger than that in the harder 50-100 keV band {SAtl'^^ x 10"^ 
erg cm~^). This, and the low peak energy, makes XRF 050416A an XRF. A detailed analysis of 
the gamma-ray properties of XRF 050416A is presented in Sakamoto et al. (2006). 

An X-ray afterglow was identifled by Kennea et al. (2005) and is discussed in detail in Mangano 



et al. (2006). The X-ray light curve follows the canonical X-ray light curve described by Nousek 
et al. (2006) with an early-time decay of ax = —2.4 it 0.5 out to 172 s, where the flux density, 
fy, is related to the time since the BAT trigger by fu{t) oc t". The canonical light curve has 
an initially steep decay due to curvature effects in the relativistic fireball. This is followed by a 
slow decay due to energy injection. The light curve makes a second transition to a decay slope 
of approximately —1 when energy injection ends. The slow decay section for XRF 050416A has 
a slope of ax = —0.44 it 0.13 and lasts from 172 to 1450 s. The late-time decay has a slope of 
ax = —0.88 lb 0.02, which continues until at least 42 days after the BAT trigger. The X-ray 
spectrum after 1450 s has a slope of (5x = — 1-04]']q5|, where the relationship between flux density 
and frequency is fy oc v^ , and shows evidence for a hydrogen column density in the host galaxy 
of Nfj = 6.8l]^'2 X 10^^ cm^^. If we assume the Predehl & Schmitt (1995) relationship between 
hydrogen column density and extinction in the Milky Way, Xjj = (1.79 x \(P''^)Av-, then this column 
density implies Ay = 'iS'^^j in the host galaxy along the line of sight to XRF 050416A. 

The optical afterglow was first identified by Cenko & Fox (2005) using the Palomar 200-inch 
Hale Telescope and confirmed by Anderson et al. (2005) using the ANU 2.3-m telescope. Fox 
(2005) noted that the afterglow is visible in the UVOT's UVW2 filter and stated that this indicates 
a redshift of z < 1. A spectrum of the host galaxy was obtained 51 days after the burst by Cenko 
et al. (2005). They measured a redshift for the host of z = 0.6535 ±0.0002 from [Oil], H/?, H7, and 
H(5 emission lines. The host's blue colour suggests ongoing star formation, as is seen in many GRB 
host galaxies. UVOT magnitudes, based on preliminary photometric calibrations, were published 
by Schady et al. (2005a,b). This paper supersedes those results. Soderberg et al. (2005) found a 
4.86 GHz radio afterglow with a flux density of 260 ± 55 /t/Jy approximately 5.6 days after the BAT 
trigger. 

Excluding the short-hard bursts (Kouveliotou et al. 1993) XRF 050416A, with a redshift of 
0.6535, is one of the nearest of the Swift-detected GRBs with spectral redshift determinations. Its 
Tgo duration and Ep are similar to those of the XRF 050406 (Krimm et al. 2005; Schady et al. 2006) 
and are fairly typical of XRFs in general. Sakamoto et al. (2006) have shown that XRF 050416A is 
consistent with the relationship between the isotropic energy equivalent, -Eiso and the peak energy 
in the rest frame, E'J-, found by Amati et al. (2002) for classical long-soft GRBs. This suggests that 
XRFs and GRBs are intimately related. XRF 050416A does not, however, follow the the empirical 
relation between E' E[so, and the jet break time in the rest frame, t', of Liang & Zhang (2005), 
which suggests that either this XRF has an unusually late jet break or that something is preventing 
us from seeing the jet break, such as a long-lasting emission component. Something similar was seen 
in GRB 021004, where identifying the jet break was complicated by several emission components 
(Holland et al. 2003). 

In this paper we present space- and ground-based ultraviolet, optical, and infrared observations 
of XRF 050416A. We have adopted a cosmology with a Hubble parameter of Ho = 70 km s~^ Mpc^^, 



^Throughout this paper we use primes to indicate quantities in the rest frame. 



a matter density of 0,m = 0.3, and a cosmological constant of J7a = 0.7. For this cosmology a 
redshift of z = 0.6535 corresponds to a luminosity distance of 3917 Mpc and a distance modulus of 
42.96. One arcsecond corresponds to 11.48 comoving kpc, or 6.95 proper kpc. The look back time 
is 6.05 Gyr. 



2. Observations 

The optical afterglow for XRF 050416A is located at R.A. = 12:33:54.6, Dec. = +21:03:27 
(J2000) (Cenko & Fox 2005), which corresponds to Galactic coordinates of (6", /") = (+82?7316, 268?7316). 
The reddening maps of Schlegel et al. (1998) give a Galactic reddening of Eb-v = 0.03 ± 0.02 mag 
in this direction. The corresponding Galactic extinctions are Au = 0.16, Ab = 0.13, Ay = 0.10, 
^Rc ~ 0.08, Ai^ = 0.06, and Aks = 0.01. The Galactic extinctions in the UVOT ultraviolet filters 
were calculated using the Milky Way extinction law from Pei (1992). The ultraviolet extinctions 
are ^uvwi = 0.23, ^UVM2 = 0.29, and AuvW2 = 0.21. 

The field of XRF 050416A is shown in Figure 1. 



2.1. UVOT Data 

The Swift spacecraft slewed promptly when the BAT detected XRF 050416A, and UVOT 
began imaging the field 65 s after the BAT trigger. A range of 14-16 exposures were taken in 
each UVOT fiher between 11:05:49 UTC on 16 Apr 2005 and 03:45:46 UTC on 18 Apr 2005. 
The UVOT has UBV filters which approximate the Johnson system, and three ultraviolet filters: 
UVWl with a central wavelength of Aq = 2600 A, UVM2 with Aq = 2200 A, and UVW2 with 
Ao = 1930 A. Examination of the first settled observation (a 100 s full-frame exposure in V) 
revealed a new source relative to the Digital Sky Survey inside the XRT error circle. This source 
had a magnitude oi V = 19.19 ± 0.29. The UVOT position of the source is R.A. = 12:33:54.596, 
Dec. = +21:03:26.07 (J2000) with a statistical error of O'.'Ol and an absolute astrometric accuracy 
of 0'.'56 (90% containment). This position is 0'.'3 from the reported XRT position (Mangano et al. 
2006). It is inside the XRT error circle and consistent with the ground-based detection reported 
by (Cenko &: Fox 2005). Subsequent exposures showed the source to be rapidly fading. 

The UVOT took 342 exposures of the field containing XRF 050416A between 2005 Apr 16 
and 2005 May 13 UTC. For the vast majority of these exposures the afterglow was too faint to be 
detected. The UVOT photometry that was used in this paper is shown in Figures 2, 3, and 4 and 
presented in Table 1. Magnitudes have not been corrected for extinction. All upper limits reported 
in Table 1 are 3a and were obtained using a 2" radius aperture with aperture corrections applied 
as discussed below. 
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Fig. 1. — This figure shows the 2 x 1200 second DFOSC y-band image taken 60 415 seconds after 
the BAT trigger. The large circle is the 3' radius BAT error circle (Sakamoto et al. 2005b). The 
small circle is the 3'.'3 radius XRT error circle (Mangano et al. 2006). North is up and east is to 
the left. 



-a 

-i-j 
.1 — I 





1 


1 1 II 1 II 


Mill 


1 1 1 1 1 1 II 


1 1 1 1 1 1 II 


1 1 1 1 1 1 II 1 1 1 1 1 1 II 


18 


— 


_ 


_ — 






GRB 050416A Jl 


20 


— 






— ~ 




— 


22 


— 


UVW2 






A 




— 


24 








w 




- 


18 


— 














— 


20 


— 










— 


22 
24 


1 


1 1 1 1 


/M2 

1 1 1 






f- 


1 1 1 II 1 1 1 


1 1 1 1 1 1 1 1 


18 


1 


1 1 1 1 


1 1 1 


— 






1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 


20 


— 






— — 




— 


22 


— 


UVWl 






~ — — ~— n 


1 1 1 1 


24 


1 


1 1 II 1 II 


Mill 


1 1 1 1 1 1 II 


1 1 1 1 1 1 II 


1 1 1 II 1 II 1 1 1 1 1 1 II 



100 1000 104 105 IQ6 10'^ 

Seconds since 2005 Apr 16.4616262 UT 



Fig. 2. — This Figure shows the UVW2, UVM2, and UVWl photometry of the afterglow of 
XRF 050416A. Red circles indicates UVOT detections. Upper limits are represented by horizontal 
lines. The horizontal error bars on the detections indicate the total period during which data was 
collected, not the exposure time. 
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Fig. 3. — This Figure shows the UVOT (red, circles) and DFOSC (green, squares) U, B, and V 
photometry of the afterglow of XRF 050416A. The details are the same as for Figure 2. 
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Fig. 4. — This Figure shows the DFOSC (green, squares) and NOTCam (blue, triangles) Re, Ic-, 
and Kg photometry of the afterglow of XRF 050416A. The details are the same as for Figure 2. 



Table 1. Photometry of XRF 050416A. The time tmid is the time from the BAT trigger to the 

middle of the observation. The exposure time is denoted by At. Upper limits are 3a upper limits. 

Errors are la statistical errors and do not include the systematic errors in the photometric zero 

points. 



t„,id (s) At (s) 



Filter 



Mag 



Err 



Instrument 



t^id (s) At (s) 



Filter 



Mag 



Err 



Instrument 



126 


100 


V 


19.19 


228 


10 


B 


> 19.49 


257 


10 


V 


> 19.39 


285 


10 


UVWl 


> 18.36 


313 


10 


B 


> 19.17 


342 


10 


V 


> 18.36 


355 


10 


UVM2 


> 18.57 


369 


10 


UVWl 


> 18.61 


397 


10 


B 


> 19.44 


482 


10 


B 


> 19.39 


497 


10 


UVW2 


> 18.41 


510 


10 


V 


> 18.61 


566 


10 


B 


> 19.32 


594 


10 X 10 


U 


19.75 


650 


10 


B 


> 19.38 


666 


10 


UVW2 


> 18.20 


679 


10 


V 


> 18.13 


707 


10 


UVWl 


> 19.30 


735 


10 


B 


> 19.49 


792 


10 


UVWl 


> 18.16 


819 


10 


B 


> 19.48 


848 


10 


V 


> 18.33 


876 


10 


UVWl 


> 18.09 


904 


10 


B 


> 19.57 


919 


10 


UVW2 


> 19.05 


933 


10 


V 


> 18.22 


1035 


100 


B 


20.74 


1140 


100 


UVW2 


> 19.81 


1245 


100 


V 


> 19.55 


1349 


100 


UVM2 


> 20.75 


1454 


100 


UVWl 


> 20.32 


1558 


100 


U 


> 20.44 


1663 


100 


B 


> 20.90 


1768 


100 


UVW2 


> 19.55 


1868 


93 


V 


19.51 


2081 


100 


UVWl 


> 20.43 


2185 


100 


U 


> 20.30 


2290 


100 


B 


> 20.65 


6850 


900 


V 


> 21.02 


11730 


900 


UVWl 


> 21.81 


12 637 


900 


u 


> 21.61 


13 508 


826 


B 


> 21.71 


18 532 


900 


V 


> 20.88 


23 601 


900 


UVWl 


> 21.93 


24 508 


900 


u 


> 21.55 


25 231 


530 


B 


> 21.60 



0.29 



0.30 



0.42 



40 



UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT'' 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 

UVOT 



59 659 


900 


U 


> 21.57 






UVOT 


60166 


99 


B 


> 20.72 




UVOT 


60 415 


2 X 1200 


V 


23.13 


0.17 


DFOSC' 


63 859 


900 


B 


> 21.86 




UVOT 


67559 


9202 


UVM2 


23.38 


0.42 


UVOT'' 


68 693 


3 x600 


Ic 


21.02 


0.17 


DFOSC 


70 541 


900 


UVWl 


> 21.98 






UVOT 


70 549 


1500 


B 


23.11 





25 


DFOSC 


71393 


791 


U 


> 21.55 






UVOT 


73 345 


10146 


UVW2 


23.25 





33 


UVOT'' 


76 988 


459 


V 


> 20.53 






UVOT 


82 070 


900 


UVWl 


> 21.92 






UVOT 


82 944 


835 


U 


> 21.44 






UVOT 


86 973 


900 


B 


> 21.96 






UVOT 


88 743 


810 


V 


> 20.83 






UVOT 


93 637 


900 


UVWl 


> 22.15 






UVOT 


94 514 


841 


u 


> 21.47 






UVOT 


98 641 


900 


B 


> 22.01 






UVOT 


100162 


311 


V 


> 20.45 






UVOT 


105 335 


900 


UVWl 


> 21.96 






UVOT 


110415 


900 


B 


> 21.95 






UVOT 


112 038 


514 


V 


> 20.47 






UVOT 


117 400 


900 


UVWl 


> 22.04 






UVOT 


117 970 


223 


u 


> 20.80 






UVOT 


122 543 


900 


u 


> 21.73 






UVOT 


123 434 


868 


B 


> 21.80 






UVOT 


129 343 


624 


V 


> 20.72 






UVOT 


135 246 


390 


UVWl 


> 21.34 






UVOT 


138 805 


401 


UVWl 


> 21.45 






UVOT 


140 587 


900 


u 


> 21.60 






UVOT 


141 137 


184 


B 


> 20.91 






UVOT 


142 008 


25 X 400 


Re 


22.84 


0.11 


DFOSC 


144 815 


900 


B 


22.19 


0.47 


UVOT 


146 738 


553 


V 


> 20.55 




UVOT 


231 672 


53 X 300 


Ic 


> 23.10 




DFOSC 


313 726 


29 X 300 


Ic 


22.74 


0.19 


DFOSC 


406 438 


21 X 200 


Ic 


23.03 


0.23 


DFOSC 


868 505 


20 957 


V 


> 22.56 






UVOT" 


1039340 


8274 


V 


> 22.06 






UVOT'' 


1384870 


6755 


V 


> 21.91 






UVOT" 


1468840 


4722 


V 


> 21.69 






UVOT" 


1941288 


14167 


V 


> 22.38 






UVOT" 


2338167 


21133 


V 


> 22.58 






UVOT" 


2424904 


26 389 


V 


> 22.63 






UVOT" 


i 146 746 


39 X 52 


Ks 


> 19.57 






NOTCam" 


5 371121 


12 866 


V 


> 22.15 






UVOT" 
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Table 1 — Continued 



tmid (s) 


At(s) 


Filter 


Mag 


Err 


Instrument 


tmid (s) 


At(s) 


Filter 


Mag 


Err 


Instrument 


30 599 


900 


V 


> 20.88 




UVOT 


3 455 561 


4722 


V 


> 21.81 




UVOT» 


36 640 


858 


UVWl 


> 22.00 




UVOT 


3 590 733 


15 986 




> 22.34 




UVOT'' 


41791 


900 


U 


> 21.51 




UVOT 


3 732 868 


18 806 


y 


22.59 


0.42 


UVOT'' 


42 552 


606 


B 


> 21.61 




UVOT 


4170 209 


321 


Y 


> 19.97 




UVOT'' 


48 456 


373 


V 


> 20.44 




UVOT 


5 706 808 


25 285 




> 22.57 




UVOT" 


51023 


36 X 78 


Ks 


19.07 


0.17 


NOTCam^' 


5 793 106 


20 898 




> 22.43 




UVOT" 


52 031 


271 


V 


> 20.35 




UVOT 


5 880 015 


23 533 


y 


> 22.52 




UVOT" 


54 363 


132 


UVWl 


> 20.78 




UVOT 


6 136 666 


15 285 




> 22.21 




UVOT" 


57 086 


6x 900 


Re 


22.29 


0.08 


DFOSC' 


6 310 565 


13 933 




> 22.13 




UVOT" 


57 990 


773 


UVWl 


> 21.70 




UVOT 


6 394 649 


10 004 


\/ 


> 21.95 




UVOT" 



"This data point consists of coadded data. 
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We performed photometry on each UVOT exposure using a circular aperture with a radius 
of 2" centred on the position of the optical afterglow. This radius is approximately equal to the 
full-width at half-maximum (FWHM) of the UVOT point-spread function (PSF). The PSF varies 
with filter and with the temperature of the telescope, so we did not match the extraction aperture 
to the PSF for each exposure. The PSF FWHM, averaged over the temperature variations, ranges 
from 1'.'79±0'.'05 for the V filter to 2'.' 17 ± (/.'OS for the UVW2 filter. The background was measured 
in a sky annulus of inner radius 17'.'5 and width 5" centred on the afterglow. The source was not 
detected in the individual ultraviolet filter exposures, so these exposures were coadded. The source 
was detected in the coadded UVW2 and UVM2 exposures, but not in the coadded UVWl exposure. 

Aperture corrections were computed for each exposure to convert the 2" photometry to the 
standard aperture radii used to define UVOT's photometric zero points (6" for UBV and 12" for 
the ultraviolet filters). Six isolated stars were used to compute the aperture correction for each 
exposure. The RMS scatter in the mean aperture correction for a single exposure was typically 
~ 0.02 mag. The RMS scatter for each exposure was added in quadrature to the statistical error in 
the 2" magnitude to obtain the total la error in each point. All detections above the 2a significance 
level are tabulated in Table 1. 

The UVOT is a photon-counting device with a frame read-out time of 11.0329 ms. It is only 
able to record one photon per detector cell during each read out. This results in coincidence losses at 
high count rates. For very high count rates, corresponding to F < 13.5, these losses are significant 
and can dramatically affect the photometry, so coincidence loss corrections must be made. We have 
corrected all of our data for coincidence loss, however the afterglow has V > 19, so coincidence 
losses are negligible, typically less than 0.01 mag, for the afterglow. Coincidence loss corrections, 
however, are significant for the stars used to compute the aperture corrections. 

The zero points used to transform the instrumental UVOT magnitudes to Vega magnitudes 
were taken from in-orbit measurements as obtained from the HEASARC Swift/UVOT Calibration 
Database (CalDB)^ dated 2005-11-18 in the file swuphot20041120v102.fits. Colour terms were 
not applied to the photometric calibrations, but preliminary calibrations of on-orbit data suggest 
that they are negligible. 



2.2. Ground-Based Data 

2.2.1. DFOSC 

The XRF 050416A optical afterglow was observed with the 1.54m Danish telescope equipped 
with DFOSC. DFOSC is focal reducer camera based on a backside illuminated EEV/MAT CCD, 
providing a pixel scale of 0'.'39 per pixel and a field of view of 13f6 x 13'6. The DFOSC observations 



http://swift.gsfc.nasa.gov/docs/hoasarc/caldb/swift/ 
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were carried out on five consecutive nights starting on 2005 Apr 17.1 UTC (see Table 1 for further 
details). Data reduction was performed following standard procedures (bias level subtraction and 
sky flat fielding) running under IRAF^. 

Photometry was performed on the combined DFOSC images using the SExtractor software 
(Bertin & Arnouts 1996) with its default parameters. Magnitudes and magnitude errors were 
determined using the MAG_AUTO option. The DFOSC data were calibrated using the BVRcIc 
field photometry of Henden (2005). All of the stars that were in common between the Henden (2005) 
photometry and our SExTRACTOR photometry were matched and used to compute photometric 
zero points for each DFOSC image. Colour terms were not used as they did not improve the quality 
of the calibration. 



2.2.2. NOT Cam 

Near-infrared observations of XRF 050416A were taken on 2005 Apr 17 and 2005 May 22 with 
NOTCam mounted on the Nordic Optical Telescope. NOTCam is a multi-mode instrument that 
is based on a 1024 x 1024 HgCdTe "HAWAII" detector providing a scale of 0'.'233 per pixel and a 
field of view of 4' x 4'. Data reduction was done by following a standard procedure for near- infrared 
imaging (shift and expose, subtract a scaled median sky, divide by twilight flat, align and combine 
by adaptive sigma clipping). The individual images were aligned from the centroid of the star 
12335122+2104142 in the 2MASS catalogue (Skrutskie et al. 2006). 

Photometry was performed on the NOT images using SExTRACTOR with the default param- 
eters in the same way as for the DFOSC images. The NOTCam photometry was calibrated using 
star 12335122+2104142. The 2MASS catalogue magnitude of this star is Ks = 11.383 ± 0.018 
(Skrutskie et al. 2006). Our infrared photometry is listed in Table 1. 



3. Results 

3.1. The Decay Rate 

We use our photometry to constrain the rate of a power law decay of the afterglow in each filter 
where there is multi-epoch photometry. The X-ray data (Mangano et al. 2006) shows a change in 
the X-ray light curve at 1450 seconds which Mangano et al. (2006) interpret as the end of energy 
injection. In order to directly compare the ultraviolet/optical/infrared decay rate to the X-ray 
decay rate we only consider photometry with mid-point times more than 1450 s after the BAT 
trigger. This is the period that Mangano et al. (2006) call Phase C and corresponds to the time 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities with National Science Foundation. 
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after the second X-ray break. 

We find UB = +1.18 ± 0.68, ay = -0.96 ± 0.53, aR = -0.57 ± 0.14, and a/ = -1.04 ± 0.66. 
It is not clear if the increase in the B-hanA flux at 144 815 s (~ 1.7 days) is real or a noise spike 
causing the source to appear brighter than it really is. The mean decay rate, excluding the i?-band 
data, is a = —0.86 ± 0.15 (standard error). This is in good agreement with the X-ray decay rate 
of ax = —0.88 lb 0.02 at the same time. Further, the decay rates in each filter, excluding i?, are 
all within 2o" of the X-ray value. This suggests that, for t > 1450 s, there is no spectral break 
between the X-ray and optical bands. As a further test we used the spectral energy distribution 
from Section 3.2 to convert all the photometry to the Re band and fit a single power law. This 
yielded a slope of —0.75 it 0.19, in good agreement with the mean slope that we determined above. 
We prefer to adopt the mean slope since it does not depend on uncertainties, and possible temporal 
variations, in the spectral shape. 

The early time UVOT y-band photometry is consistent with the afterglow having a constant 
F-band magnitude between 126 and 1868 seconds after the BAT trigger. The decay rate in this 
interval is ay = — O.llibO.17. This is consistent with the X-ray slope during Phase B (—0.44+0.13) 
(Mangano et al. 2006) at the 1.6cr level. 



3.2. The Spectral Energy Distribution 

We constructed the spectral energy distribution between 1930 A and 22 200 A at 60 415 s after 
the BAT trigger. The photometry data nearest this time for each filter were transformed to 60 415 
s assuming a decay rate of a = —0.88. We adopted the X-ray decay rate instead of the optical 
decay rate that we determined in Section 3.1 because the X-ray decay is better constrained, and 
we believe that the X-ray and optical decays are the same (see Sect. 3.1). The UVW2 and UVM2 
magnitudes were transformed to flux densities using the conversion factors in the Swift/UVOT 
CalDB. The optical and infrared magnitudes were converted to flux densities using the zero points 
of Fukugita et al. (1995) and Cox (2000). Each data point was corrected for Galactic extinction 
using the reddening value of Schlegel et al. (1998) (see Section 2) but not for any extinction that 
may be present in the host galaxy or in intergalactic space along the line of sight to the burst. The 
SED is shown in Figure 5. 

The SED was flt by fu{v) oc i/f^lO^^^^^'^' where fu{v) is the flux density at frequency v, f5 is 
the intrinsic spectral index, and A{v) is the extragalactic extinction in the host {z = 0.6535) at 
frequency v. We found that the best fit occurred for the Small Magellanic Cloud (SMC) extinction 
law of Pei (1992), although using a Milky Way or Large Magellanic Cloud extinction law does not 
significantly change the results. We prefer the SMC since it provides good fits to the extinction seen 
in the host galaxies of other GRBs (e.g., Holland et al. 2003). If we fix the intrinsic spectral slope to 
the X-ray value of [3x = —1.04 (Mangano et al. 2006) the extinction in the host is Ay = 0.24 + 0.06 
mag. If we allow both (3 and Ay to be free parameters the best fit occurs for [3 = — 1.14 + 0.20 and 
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Fig. 5. — This Figure shows the SED of the ultraviolet/optical/infrared afterglow of XRF 050416A 
at 60 415 s (0.7 days) after the BAT trigger. The filled circles represent observed photometry 
corrected for extinction in the Milky Way. The solid line represents the SED fit with an extinction 
in the host of Ay = 0.19 it 0.11 mag assuming an intrinsic power law spectrum with a slope of 
j3 = —1.14 zb 0.20. The dashed line shows the unreddened spectrum. 
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Ay = 0.19 ± 0.11. Therefore, we believe that there is approximately 0.2 mag of F-band extinction 
in the host along the line of sight to XRF 050416A. 

The low value of the extinction in the host is at odds with the high neutral hydrogen column 
density found by fitting the X-ray spectrum (Mangano et al. 2006). Predehl & Schmitt (1995) find 
A^H = (1.79 X 10^^)yly for the the conversion between hydrogen column density and extinction in 
the Milky Way. Using the value of A^h derived from the Swift/XRT X-ray observations (Mangano 
et al. 2006) we find Ay = S.SJIq'^, which implies significant extinction in the host along the line 
of sight to XRF 050416A. The discrepancy between this and the small amount of extinction that 
we find from the ultraviolet/optical/infrared SED suggests that the relationship between hydrogen 
gas and dust is different in the host of XRF 050416A than it is in the Milky Way. Combining the 
A-ray hydrogen column density with our estimate of Ay gives Ah = (36_22 x 10^^)yly in the host 
galaxy. This is consistent with the Nn/Ay relation found in the SMC of Ajj = (15.4 x lO'^^) Ay 
using Eq. 4 and Table 2 of Pei (1992). 

The high gas-to-dust ratio seen in the host of XRF 050416A is typical of other GRB host galax- 
ies, such as GRB 000301C (Jensen et al. 2001), GRB 000926 (Fynbo et al. 2001), and GRB 020124 
(Hjorth et al. 2003), all of which have ratios consistent with that observed in the SMC, (see also 
Stratta et al. 2004; Kann et al. 2006). This suggests that our choice of using an SMC extinction 
law to determine the intrinsic spectral slope is reasonable. The high gas-to-dust ratio may be a 
consequence of dust destruction by the ultraviolet and A-ray fiux from the GRB as described by 
(e.g. Waxman & Draine 2000; Perna et al. 2003). Alternately, it may indicate that the star forma- 
tion in the host is fairly recent and that there may not have been time for large amounts of dust 
to form. 



3.3. Cooling Break 

The VRqIc photometry yields a mean decay index of q = —0.86 it 0.15 for times more than 
1450 seconds after the burst. The ultraviolet /optical/infrared SED at 60 415 s has a slope of 
(3 = —1.14 lb 0.20, after correcting for extinction in the host galaxy. This is consistent with the 
intrinsic spectral slope derived from the A-ray data. The lack of evidence for a change in either 
the intrinsic spectral slope or the decay rate between the A-ray and optical bands suggests that 
there is no spectral break between these two regimes. Therefore the synchrotron cooling frequency 
must lie either above the A-ray band or below the Ks band at 60 415 s (~ 0.7 days) after the BAT 
trigger. 

The spectral and temporal decay rates can be used to predict the index of the electron energy 
distribution, p. The decay rate will be a more robust number than the spectral index since the 
spectral slope may be affected by uncertainties in the absorption along the line of sight. After the 
jet break the magnitude of the decay index should equal the value of the electron index. However, 
electron indices of less than unity are unphysical. Therefore, we can rule out a jet break before 42 
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days after the BAT trigger because the X-ray decay rate remains constant at —0.88 between 1450 
s and at least 42 days. 

In order to determine the location of the cooling break, and to estimate the value of the electron 
index, we use closure relations between oi and /?. These are based on the relations of Dai & Cheng 
(2001) when 1 < p < 2. For p > 2 the relations of Sari et al. (1999) were used for the case of a 
homogeneous circumburst environment, and the relations of Chevalier & Li (1999) were used for 
the case of a pre-existing stellar wind. Table 2 lists the closure relations and their observed values. 
None of the closure relations for p > 2 give values consistent with zero. Therefore, we believe that 
1 < p < 2. This is a low value for a GRB, but three of the ten bursts studied by Panaitescu & 
Kumar (2002) had p ~ 1.4, so it can not be ruled out. For 1 < p < 2 the case of a homogeneous 
circumburst medium with vx < i^c, and the cases of a pre-existing stellar wind with i^x < t'c or 
J^c < t'opt ah give comparable closure values that are consistent with zero. 

If 1 < p < 2 then the predicted intrinsic spectral slope is between and —0.5 if the cooling 
break is above the X-ray band and between —0.5 and —1 if the cooling break is below the optical 
bands. Both the X-ray and ultraviolet/optical/infrared spectral slopes are consistent with /3 = —1, 
which corresponds to p = 3 for i^x < t'o which is unusually high for a GRB. However, if i/c < t'opt 
then p = 2, which is fairly typical for GRBs. An electron index of p = 2 is consistent with cases 5 
and 7 of the closure relations in Table 2. In these two cases the predicted optical decay is a = — 1 
regardless of the density structure of the environment that the burst is expanding into. This decay 
slope is slightly steeper than the observed decay rate of —0.86 ± 0.15, but it is within the observed 
uncertainties. Therefore, we believe that p ~ 2 and Vc < '^opt for XRF 050416A at 60 415 s after 
the BAT trigger. 



3.4. Energy Considerations 

Sakamoto et al. (2006) find that the isotropic equivalent energy of XRF 050416A is 1.2 x 10^^ 

Table 2. This Table lists the closure relationships for the various cases under consideration. The 

theoretical value of each closure relation is zero. 



Case 


Model 


Environment 


P 




Closure 


Value 


1 


I'opt < l^X < I^c 


ISM 


l<p< 2 


a - 


-3/8/3 + 9/16 


+0.13 ±0.19 


2 






p>2 


a - 


-3/2/3 


+0.85 + 0.11 


3 




Wind 


1<P<2 


a - 


- 1/4/3 + 9/8 


+0.55 + 0.16 


4 






p>2 


a - 


- 3/2/3 + 1/2 


+0.35 + 0.11 


5 


l^c < I'opt < I^X 


ISM 


l<p<2 


a - 


- 3/8/3 + 5/8 


+0.19 + 0.19 


6 






p>2 


a - 


- 3/2/3 - 1/2 


+0.55 + 0.11 


7 




Wind 


l<p< 2 


a - 


- 1/4/3 + 3/4 


+0.18 + 0.16 


8 






p>2 


a - 


- 3/2/3 - 1/2 


+0.55 + 0.11 
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erg. This, and the observed hmits on the jet break time, can be used to estimate the opening angle 
of the jet and thus the total gamma-ray energy of the burst (Rhoads 1999; Sari et al. 1999; Frail 
et al. 2001). The jet opening angle for XRF 050416A is 

% = O.lOetf ^(??^/0.2)^/^(n/0.1)^/^ (1) 

where tj is the observed jet break time in days since the BAT trigger, r/^ is the efficiency of converting 
energy in the ejecta into gamma rays, and n is the particle density in cm^'^. The X-ray light curve 
(Mangano et al. 2006) suggests that there is no jet break at times earlier than 42 days. Setting 
tj > 42 yields 9j > 25° assuming rj^ = 0.2 and n = 0.1 cm~^. This result is not strongly sensitive 
to our choices of rj^y and n. The corresponding energy in gamma rays, corrected for beaming, is 
E'y > 1.1 X 10^*^ erg. The upper limit on the gamma-ray energy can be obtained from the case 
where there is no beaming (and thus no jet break). In that case E^ = E^gQ, which implies that 
1.1 X 10^'^ erg < Ey < 1.2 x 10^^ erg. Therefore the total energy must be less than the canonical 
GRB energy of Bloom et al. (2003) if there is any beaming. 

Sakamoto et al. (2006) showed that XRF 050416A fits onto the Amati relation (Amati et al. 
2002) between E^g^ and E' . They also showed that this burst can not be made to satisfy the 
Ghirlanda relation (Ghirlanda et al. 2004) between E^ and E' due to the lack of a jet break out 
to 35 days after the BAT trigger. Mangano et al. (2006) have shown that there is no jet break to 
at least 42 days, which makes the discrepancy with the Ghirlanda relation even larger. 

The Liang & Zhang (2005) relation 

0-85 x—£- X -L- (2) 



1053 erg ■ VlOOkeVy V^day 

implies that the jet break should occur at ~ 1.5 days after the BAT trigger in the observer's frame. 
If we assume that XRF 050416A has a two-component jet with a narrow component that breaks 
at 1.5 days the observed decay between 1.5 and 42 days will be the sum of the post-narrow-jet 
break component for the narrow jet and the post-cooling break component for the wide jet. A 
two-component jet model has been proposed to explain some of the observed features in the decays 
of GRB 030329 (Berger et al. 2003; Sheth et al. 2003), such as the presence of two achromatic 
breaks. We propose that the observed optical afterglow of XRF 050416A may be the sum of two 
jets that both conform with the Ghirlanda relation and the Liang & Zhang relation. In this picture 
the narrow jet would have an opening angle of Oj^n ~ 7° assuming tj^n = 1-5 days, and the wide jet 
would have 6j^w ^ 25° for a wide jet break time of tj^w > 42 days. The gamma-ray energies in each 
jet are E^^n ~ 5 x 10^^ erg and E^^yj ^ 7 x 10^^ erg, which is roughly consistent with the division 
of energy in the two-component jet model for GRB 030329. 

Peng et al. (2005) find that for two component jets where the total energy in the wide jet is 
greater than that in the narrow jet the wide jet will dominate the optical afterglow after approx- 
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imately 0.1-1 days. Since the deceleration time of the wide jet is similar to the break time of the 
narrow jet the emergence of the wide jet at this time can mask the steepening of the light curve 
caused by the jet break in the narrow component. We suggest that this is what has happened in 
XRF 050416A. We see no evidence for a change in the spectral slope or the decay rate between the 
optical and X-ray regimes at 60415 s (~ 0.7 days). This implies that both the narrow and wide 
components have cooling frequencies below the optical band at this time. 



3.5. Constraints on a Supernova Component 

UVOT observation in the y-band filter were taken up to 74 days after the BAT trigger. We 
searched the late time exposure for evidence of a rebrightening. There is one low-significance UVOT 
detection of the afterglow after 144 815 s (1.676 d). There is a second low-significance detection at 
3 732 868 s (43 d) . A visual examination of the second detection suggests that it is a noise spike in 
the data and not a point source. 

A type Ib/c supernova like SN1998bw (Patat k Piemonte 1998) at the distance of XRF 050416A 
is expected to peak at approximately 10(1 + z) ^ 16 days after the burst. We find no evidence for a 
source at the location of the afterglow at this time down to a 3a limiting magnitude of Vjim = 21.9 
in coadded exposures. For z = 0.6535 the observed V band approximately corresponds to the rest 
frame U band, so our upper limit corresponds to Mjj < —21.1. SN1998bw had a peak t/-band ab- 
solute magnitude of Mjj = —19.16 (Galama et al. 1998) so the UVOT data is not able to constrain 
the existence of a supernova component to the afterglow of XRF 050416A. 



4. Conclusions 

XRF 050416A is a near-by {z = 0.6535), short (Tgo = 2.4 s) XRF that shows no evidence for a 
jet break out to at least 42 days after the burst. The spectral slope and decay rate are the same in 
the optical as in Phase C of the X-ray decay (Mangano et al. 2006). This suggests that there is no 
coohng break between approximately 1.24 A (10 keV) and 22 000 A (0.006 keV) at 60 415 s after 
the burst. Further, the constancy of the decay suggests that the cooling break did not pass through 
the optical between 1450 s (0.016 days) and 406 438 s (4.7 days), nor did it pass through the X-ray 
regime between 1450 s (0.016 days) and ~ 3.6 x 10^ s (42 days). We find the best agreement with 
the synchrotron model occurs if the cooling break is below the optical at 60 415 s (0.7 days) after 
the burst. 

The optical light decay slope is a = —0.86 ± 0.15 and the intrinsic ultraviolet/optical/infrared 
spectral slope is /? = —1.14 it 0.20. The best estimate of the electron index is p ~ 2. We are 
unable to distinguish between a burst occurring in a homogeneous environment or a wind-stratified 
one. The lack of a jet break out to > 42 days implies that the jet opening angle is > 25° and the 
gamma-ray energy is 1.1 x 10^'^ erg < E^ < 1.2 x 10^^ erg. If XRF 050416A has the canonical 
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energy of Bloom et al. (2003) then the burst is not beamed. 

An alternate interpretation of the data is that XRF 050416A has a two component jet. The 
narrow component has an opening angle of 9j^n ~ 7° and experienced a jet break at 1.5 days, in 
agreement with the prediction of the Liang & Zhang (2005) relation between £^iso) E' and t'-. The 
wide component does not break until at least 42 days after the burst and must have an opening 
angle of > 25°. The wide component jet contains approximately 14 times more energy than the 
narrow component jet does. This behaviour is similar to what was seen in GRB 030329 (Berger 
et al. 2003; Sheth et al. 2003). 

We find no evidence for a large amount of extinction along the line of sight to XRF 050416A. 
The best fit ultraviolet/optical/infrared spectrum suggests that Ay ~ 0.2 mag in the host galaxy. 
This is inconsistent with the large hydrogen column density implied by the X-ray spectrum. If we 
assume that the gas-to-dust ratio in the host is the same as that in the Milky Way the implied 
extinction is Ay = S.S'^q'j mag. Using the extinction derived from the optical data and the 
hydrogen column density derived from the X-ray data the gas-to-dust ratio in the host is N^/Ay = 
3.6^2 2 ^ 10^^. This is dramatically larger than it is in the Milky Way, but consistent with what is 
seen in the SMC, and in the host galaxies of other GRBs. 

The authors wish to thank Scott Barthelmy, and the GRB Coordinates Network (GCN) for 
rapidly providing precise GRB positions to the astronomical community. This research has made 
use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, under contract with NASA. This publication 
makes use of data products from the Two Micron All Sky Survey, which is a joint project of the 
University of Massachusetts and the Infrared Processing and Analysis Center/California Institute 
of Technology, funded by the National Aeronautics and Space Administration and the National 
Science Foundation. We acknowledge the use of public data from the Swift data archive. This 
paper is based, in part, on observations taken with the Nordic Optical Telescope, operated on the 
island of Santa Miguel de la Palma jointly by Denmark, Finland, Iceland, Norway, and Sweden in 
the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. 
STH would like to thank the DARK Cosmology Centre for its hospitality while writing part of 
this paper. The DARK Cosmology Centre is funded by the Danish National Research Founda- 
tion. This work is sponsored at Penn State University by NASA's office of Space Science through 
contract NAS5-00136 and at MSSL by funding from PPARC. The research of JG is supported by 
the Spanish Ministry of Science and Education through programmes ESP2002-04124-C03-01 and 
AYA2004-01515. 



REFERENCES 

Amati, L., Fronteral, F., Tavani, M., in't Zand, J. J. M., Antonelli, A., Costa, E., Feroci, M., 
Guidorzi, C, 2002, A&A, 390, 81 



20 



Anderson, G., Salvo, M., Rich, J., & Schmidt, B. P., 2005, GCN Circular 3266 

Barthelmy, S. D., Barbier, L. M., Cummings, J. R., Fenimore, E. E., Gehrels, N., Hullinger, D., 
Krimm, H. A., Markwardt, C. B., et al., 2005, Sp. Sci. Rev., 120, 143 

Berger, E., Kulkarni, S. R., Pooley, G., Frail, D. A., Mclntyre, V., Wark, R. M., Sari, R., Soderberg, 
A. M., et al., 2003, Nature, 126, 154 

Berlin, E., & Arnouts, S., 1996, A&AS, 117, 393 

Bloom, J. S., Frail, D. A., k Kulkarni, S., 2003, ApJ, 594, 674 

Burrows D. N., Hill, J. E., Nousek, J. A., Kennea, J. A., Wells, A., Osborne, J. R, Abbey, A. F., 
Beardmore, A., et al., 2005, Sp. Sci. Rev., 120, 165 

Cenko, S. B., & Fox, D. B., 2005, GCN Circular 3265 

Cenko, S. B., Kulkarni, S. R., Gal- Yam. A., Berger, E., 2005, GCN Circular 3542 

Chevalier, R. .A. & Li, Z.-Y., 1999, ApJ, 520, L29 

Cox, A. N., ed., 2000, Allen's Astrophysical Quantities (4th ed; New York: Springer) 

Dai, Z. G., k Cheng, K. S., 2001, ApJ, 558, L109 

Fox, D. B., 2005, GCN Circular 3408 

Frail, D. A., Kulkarni, S. R., Sari, R., Djorgovski, S. G., Bloom, J. S., Galama, T. J., Reichart, 
D. E., Berger, E., et al., 2001, ApJ, 562, L55 

Fukugita, M., Shimasaku, K., k Ichikawa, T., 1995, PASP, 107, 94 

Fynbo, J. P. U., Gorosabel, J., Dall, T. H., Hjorth, J., Pedersen, H., Andersen, M. I.; M0ller, P., 
Holland, S. T., et al., 2001, A&A, 373, 796 

Fynbo, J. P. U., Jensen, B. L., Hjorth, J., Wiersema, K., Starling, R., Vreswijk, P., Rol, E., Levan, 
A., et al., 2005b, GCN Circular 3176 

Galama, T. J., Vreeswijk, P. M., van Paradijs, J., Kouveliotou, C, Augusteijn, T., Bohnhardt, H., 
Brewer, J. P.; Doublier, V., et al., 1998, Nature, 395, 670 

Gehrels, N., Chincarini, G., Ciommi, P., Mason, K. O., Nousek, J. A., Wells, A. A., White, N. E., 
Barthelmey, S. D., et al., 2004, ApJ, 611, 1005 

Ghirlanda, G., Ghisellini, G., k Lazzati, D., 2004, ApJ, 616, 331 

Henden, A., 2005, GCN Circular 3454 



21 



Hjorth, J., M0ller, P., Gorosabel, J., Fynbo, J. P. U., Toft, S., Jaunsen, A. O., Kaas, A. A., 
Purismo, T., Torii, K., et al., 2003, ApJ, 597, 699 

Holland, S. T., Weidinger, M., Fynbo, J. P. U., Gorosabel, J., Hjorth, J., Pedersen, K., Mendez, 
A. J., Augusteijn, T., et al., 2003, AJ, 125, 2291 

Jensen, B. L., Fynbo, J. P. U., Gorosabel, J., Hjorth, J., Holland, S. T., M0ller, P., Thomsen, B., 
Bjornsson, G., et al., 2001, A&A, 370, 909 

Kann, D. A., Klose, S., k Zeh, A., 2006, ApJ, 641, 993 

Kennea, J. A., Racusin, J. L., Burrows, D. N., Mangano, V., Sakamoto, T., & Gehrels, N., 2005, 
GCN Circular 3268 

Krinim, H., Barbier, L., Barthelmy, S., Cummings, J., Fenimore, E., Gehrels, N., Hullinger, D., 
Lamb, D. Q., et al., 2005, GCN Circular 3183 

Kouveliotou, C, Meegan, C. A., Fishman, G. J., Bhat, N. P., Briggs, M. S., Koshut, T. M., 
Paciesas, W. S., & Pendleton, G. N., 1993, ApJ, 413, LlOl 

Liang, E., & Zhang, B., 2005, ApJ, in press, astro-ph/0504404i 

Mangano, V., La Parola, V., Cusumano, G., Mineo, T., Campana, S., Capalbi, M., Chincarini, G., 
Giommi, P., et al., 2006, ApJ, submitted, astro-ph/0603738 

Nousek, J. A., Kouveliotou, C, Grupe, D., Page, K., Granot, J., Ramirez-Ruiz, E., Patel, S. K., 
Burrows, D. N., et al., 2006, ApJ, in press, astro-ph/0508332 

Panaitescu, A., k Kumar, P., 2002, ApJ, 571, 779 

Patat, F., k Piemonte, A., 1998, lAUC 6918 

Pei, Y., 1992, ApJ, 395, 130 

Peng, F. Konigl, A., k Granot, J., 2005, ApJ, 626, 966 

Perna, R., Lazzati, D., k Fiore, F., 2003, ApJ, 585, 775 

Predehl, P., k Schmitt, J. H. M. M., 1995, A&A, 293, 889 

Rhoads, J. E., 1999, ApJ, 525, 737 

Roming, P. W. A., Kennedy, T. E., Mason, K. O., Nousek, J. A., Ahr, L., Bingham, R. E., Broos, 
P. S., Carter, M. J., et al., 2005, Sp. Sci. Rev., 120, 95 

Sakamoto, T., Barbier, L., Barthelmy, S. D., Burrows, D. N., Chester, M., Cummings, J., Fenimore, 
E., Gehrels, N., et al., 2005a, GCN Circular 3264 



22 



Sakamoto, T., Barbier, L., Barthelmy, S. D., Cummings, J. R, Fenimore, E., Gehrels, N., Hullinger, 
D., Krimm, H., et al., 2005b, GCN Circular 3273 

Sakamoto, T., Barbier, L. Barthelmy, S. D., Cummings, J. R., Fenimore, E. E., Gehrels, N., 
Hullinger, D., Krimm, H. A., et al., 2006, ApJ, 636, L73, 

Sari, R., Piran, T., k Halpern, J. P., 1999, ApJ, 519, L17 

Schady, P., Sakamoto, T., McGowan, K., Boyd, P., Roming, P., Nousek, J. & Gehrels, N., 2005a, 
GCN Circular 3276 

Schady, P., Sakamoto, T., McGowan, K., Boyd, P., Roming, P., Nousek, J. &: Gehrels, N., 2005b, 
GCN Circular 3280 

Schady, P., Mason, K. O., Osborne, J. P., Page, M. J., Roning, P. W. A., Still, M., Zhang, B., 
Blustin, A. J., et al., 2006, AJ, submitted, astro-ph/0601182 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M., 1998, ApJ, 500, 525 

Sheth, K., Frail, D. A., White, S., Das, M., Bertoldi, F., Walter, F., Kulkarni, S. R., & Berger, E., 
2003, ApJ, 595, L33 

Skrutskie, M. F., Cutri, R. M., Stiening, R., Weinberg, M. D., Schneider, S., Carpenter, J. M., 
Beichman, C, Capps, R., et al., 2006, AJ, 131, 1163 

Soderberg, A. M., 2005, GCN Circular 3318 

Stratta, G., Fiore, F., Antonelh, L. A., Piro, L., & De Pasquale, M., 2004, ApJ, 608, 846 

Waxman, E., & Draine, B. T., 2000, ApJ, 537, 796 



This preprint was prepared with the AAS lATpjjX macros v5.2. 



